Excluded volume effects on semiflexible ring polymers 
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Two-dimensional semiflexible polymer rings are studied both by imaging circular DNA adsorbed 
on a mica surface and by Monte Carlo simulations of phantom polymers as well as of polymers 
with finite thickness. Comparison of size and shape of the different models over the full range of 
flexibilities shows that excluded volume caused by finite thickness induces an anisotropic increase 
of the main axes of the conformations, a change of shape, accomplished by an enhanced correlation 
along the contour. 



DNA is one of the most versatile building blocks for 
the self-assembly of nanoscale structures pQ. The com- 
plementary base-pairs enables almost unlimited possibil- 
ities in designing highly tailored constructs [2] . But even 
without specially designing sequences for specific base- 
pairing, biopolymers such as DNA or cytoskeletal fila- 
ments self- assemble in their natural environment. DNA 
occurs both in linear and circular forms and condenses 
into toroidal structures [3] 0] . Actin assembles into fil- 
aments and bundles [5] and also builds ring [MS], and 
racket like [H [10] complexes. From a geometrical per- 
spective filaments denote the simplest building blocks. 
The next level of complexity is reached when topological 
constrained forms as in rings arise. Considering poly- 
mer rings as building blocks their size and shape are of 
eminent importance. Both size and shape strongly de- 
pend on two internal biopolymer properties: the ability 
to bend and the effective diameter of the polymer. In- 
deed the shape of polymer rings has been investigated 
theoretically regarding the influence of the polymers flex- 
ibility L/lp, given by the ratio of total polymer length L 
and its persistence length l p [11 . However, our experi- 
mental observations show that the previously neglected 
finite thickness not only regulates the absolute size [12] 
of a polymer configuration but also its shape. This ob- 
servation is especially important in confined geometries 
utilized in the preparation of biopolymer assemblies of 
higher order [TO] . 

Coarse grained polymer models rely on phantom 
chains, which allow segments to overlap. To describe 
real polymers with finite thickness, the excluded volume 
of a polymer chain is accurately accounted for by tube 
models [13], where the tube imposes a hard core poten- 
tial, see Fig. [T] To access the effects of finite thickness 
and topology experimentally in a well-defined setup, we 
investigate circular DNA adsorbed on a mica surface, 
which has previously been shown to obey two dimen- 
sional worm-like chain statistics p~4j [15] . Our data ver- 
ifies Flory's predicted growth in size due to polymer's 
finite thickness [12] . as has also been accomplished ex- 
perimentally in different contexts [T6HT8] and theoreti- 



cally by self-consistent and renormalization group theo- 
ries [TO] , Since DNA can be produced in different lengths, 
it serves as a model system to investigate the shape over 
the full range of flexibility, which so far has only been 
forecasted theoretically for phantom polymers in three 
dimensions [TT] , 

Here we study the effects of excluded volume caused by 
finite thickness on semiflexible polymer rings by imaging 
circular DNA on a mica surface and performing Monte 
Carlo simulations of semiflexible phantom polymers and 
polymers with finite thickness. Calibrating the polymer's 
diameter d by comparing the tangent-tangent correla- 
tion, we obtain good quantitative agreement between ex- 
periments and simulations with diameter estimates from 
polyelectrolyte theory. The finite thickness leads to an 
apparent stiffening and an increase of the principal axes 
of the polymers configurations as observed by the mean 
square radius of gyration. In contrast to Flory's pre- 
diction, the asphericity as a shape measure proves the 
growth to be anisotropic, resulting in a change of shape. 

The DNA rings without superhelicity were produced 
from nicked plasmids pSHl, pBR322, and pUC19 with 
flexibilities of L/l p = 40, 30, and 18.3, respectively. Plas- 
mid pUC19 was treated with restriction enzyme Rsal to 
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FIG. 1: Simulation snap-shot of a phantom polymer (A) 
and a polymer with finite thickness (B) on a plane. Phan- 
tom polymers bear overlaps and kink-like bending resulting 
in smaller configurations with more aspherical shapes than 
polymers with finite thickness. 
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produce three different linear fragments, and using T4 
DNA ligase, minicircles of different flexibilities were ob- 
tained L/lp = 12.4, 4.6. In order to access the trajec- 
tory of the DNA rings by Atomic Force Microscopy, the 
sample was deposited on mica, see supplementary ma- 
terial for details. Previous analysis of us [15] proofed 
the tangent-tangent correlation and the mean square ra- 
dius of gyration to obey two dimensional statistics of the 
polymer rings. 

For the Monte Carlo simulation of a semiflexible poly- 
mer rings with persistence length Z p , preceding simula- 
tion methods were customized to describe two dimen- 
sional polymers and extended to incorporate finite thick- 
ness as outlined in the supplementary material. In the 
simulations with finite thickness configurations including 
intersections of tubes of the diameter d around each seg- 
ment are rejected. Uncorrelated data sets are obtained 
by taking configurations every 10 6 Monte Carlo steps for 
phantom polymers and every 10 8 steps for simulations 
with finite thickness. Large ensembles are sampled such 
that the statistical error based on a normal distribution 
of the observable is of the size of the symbols in all figures 
shown. 

Semiflexible polymers are well described by the worm- 
like chain model, where the polymer is modeled as an 
elastic rod with bending modulus n [13 j. Represent- 
ing the polymer by a differential space curve r(s) of 
length L parametrized by an arc length s, its statis- 
tical properties are determined by the elastic energy 
H = k/2 ds [dt(s)/ds] 2 , where t(s) = dr(s)/ds is 
the tangent vector. The persistence length l p as a 
measure of the stiffness is defined by the initial decay 
of the mean tangent-tangent correlation (t(s)t(s')) = 
exp( — \s — s f \/l p ) : given by l p = for a polymer em- 
bedded in two dimensions. 

Size and shape of a polymer are comprised in the radius 
of gyration tensor, 

Qij = J dsr i (s)r j (s)-^ J dsr^s) J dsTj(s) , (1) 

whose eigenvalues Ai and A2 define the spatial extent of 
the polymer along its principal axes. The squared radius 
of gyration measures the total size of an object and hence 
is given by the sum of the two eigenvalues, 
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The criterion for the shape of a polymer is the asphericity, 
which is given by the normalized variance of Ai and A2 
[20] , yielding in two dimensions, 



(Ai - A 2 ) s 
(Ai + A 2 )2 



(3) 



The asphericity ranges between and 1; A = for the 
most spherical object in two dimensions, the ring, and 
A = 1 for the most aspherical configuration, a rod. 
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FIG. 2: Tangent-tangent correlation of a semiflexible polymer 
ring for different theoretical models and for DNA (pUC19). 
The finite thickness is calibrated to d/l p — 0.13. Incorpo- 
rating the same persistence length models with higher self- 
exclusion show enhanced correlations. 



To model real polymers, two internal parameters are 
required, the persistence length l p and the diameter d of 
the filament. For the persistence length of DNA we use 
the widely accepted value of l p = 50nm. The effective 
diameter of DNA, being a polyelectrolyte, changes in a 
predictable manner in response to its surrounding ionic 
solution, as it has been determined theoretically [21] and 
experimentally [22] . For our experimental conditions the 
ratio of diameter to persistence length is d/l p = 0.13. 
To calibrate the simulation parameters to the experi- 
mental data, the tangent- tangent correlation is an ideal 
observable reflecting the statistics along the whole con- 
tour of the polymer. The agreement between simulation 
with finite thickness and experimental data is exempli- 
fied for the tangent-tangent correlation for a flexibility 
of L/l p = 18.3 considering plasmid pUC19 in Fig. [2J 
Included in the graph are also results for a phantom 
polymer, where self-intersecting configurations are per- 
mitted, and for reasons of comparison the trivial limit 
of vanishing thickness d —> denoted self-avoiding poly- 
mer. In the latter case, only intersections of the polymer 
backbone are rejected. This limit does not describe the 
experimental data quantitatively as good as the simu- 
lations with a thickness of d/l p = 0.13 (7% increase of 
the reduced chi-square). Hence, the diameter estimate 
based on polyelectrolyte theory yields good quantitative 
agreement. 

In the tangent- tangent correlation function of Fig. [2] 
two effects are evidenced if comparison is made between 
a polymer with finite thickness and a phantom chain of 
same persistence length (solid and dashed line in Fig. [5J 
respectively), one due to the finite thickness of the poly- 
mer and the other to the circular topology. On short 
distances, the finite thickness restricts the available con- 
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FIG. 3: Squared radius of gyration (Rg) compared to the size 
of the corresponding rigid ring R 2 C — (L/2tt) 2 vs. L/l p . The 
influence of finite thickness arises in the semiflexible regime 
after an effective stiffening in accordance with analytical pre- 
dictions for stiff rings. In the flexible regime, the increase in 
size is required to model the experimental data. Error bars 
indicate the statistical errors. 



formational space, thus increasing the directional correla- 
tion. In order to respect the circular topology, the corre- 
lation function must become more negative on distances 
s ~ L/2. In summary, the polymer with finite thickness 
appears effectively stiffer. 

Similarly, this effective stiffening also affects the over- 
all size of the polymer as measured by the squared radius 
of gyration (Rg). In Fig.[3](iT^) normalized with the ra- 
dius of gyration of a rigid circle, namely (Rg) is 
plotted versus the flexibility parameter L/l p for the con- 
sidered models together with the data. Three regimes are 
discerned; in the stiff regime all simulation data for the 
radius of gyration normalized by total polymer length 
versus flexibility L/l p follow the predicted linear decay 
[23] . The stiff regime extends up to high flexibilities 
compared to open chains. For phantom polymers this 
effect can be accounted for by an effective stiffening due 
to the topological constraint of a ring [23] . For poly- 
mers with finite thickness the squared radius of gyration 
follows the analytic result for the stiff limit up to even 
higher flexibilities of approximately L/l p = 12 indicating 
a further effective stiffening due to the polymer's thick- 
ness i.e. as a result of excluded volume. This is in agree- 
ment with the observation of enhanced correlations in 
the tangent-tangent correlation. The semiflexible regime 
is a crossover region, where departing from the analyt- 
ically determined stiff limit, phantom polymers show a 
linear decay as flexibility increases suggesting an initial 
step by step excitation of higher modes. In contrast, for 
polymers with finite thickness these initial higher modes 
are suppressed, resulting in a direct transition from the 
linear decay in the stiff limit to the power law decay in 



the flexible regime. Such polymers with finite thickness 
deviate from the stiff limit to larger sizes as is also ob- 
served for three-dimensional phantom polymer rings [23 . 
Finally, in the flexible regime both models have substan- 
tially different radii of gyration. In contrast to phantom 
polymers, polymers with finite thickness show notably 
larger sizes recovering Flory's swelling effect. The scaling 
exponent agrees with Flory's predicted exponent (data 
not shown). The segments of phantom polymers overlap 
strongly to maximize entropy as flexibility permits. Pre- 
cisely those modes are, however, forbidden for polymers 
with finite thickness yielding a larger mean squared ra- 
dius of gyration. Flory's argument oversimplifies a semi- 
flexible chain of segments to an ideal gas and assumes 
that all chain segments overlap with an equal probability 
with each other. This results in a growth in size that 
is equally large along all principal axes of the polymer. 
Hence, Flory's description predicts an isotropic increase 
of the principal axes, which will be tested below when 
considering the asphericity. The experimental data in 
Fig.[3]are in agreement with the finite thickness polymers, 
providing solid evidence that a description of DNA as a 
semiflexible polymer with persistence length l p = 50 nm 
and effective diameter d/l p = 1.3 is a faithful description 
of DNA conformations. Indicated in the graph (Fig. [3| 
are only statistical errors of a Gaussian distributed ob- 
servable as a lower estimate of the statistical error. Fur- 
thermore, systematic errors may arise, first, due to the 
limited resolution of the AFM images, and second, due 
to the fact that the minicircles are not nicked and may 
thus experience a slight distortion. 

In order to test Flory's prediction, the asphericity (A) 
(Eq. |3| is plotted in Fig. [4] as a function of the flexibil- 
ity parameter L/l p and the three regimes appear again. 
Starting from a ring configuration with A = for infinite 
stiffness L/l p = 0, the asphericity grows linearly for both 
models in the stiff region due to the fact that polymers 
have an elliptical shape of increasing eccentricity, as it has 
been recently predicted by scaling arguments [IT]. Self- 
exclusion plays no role in these configurations because 
the segments are well separated from each other due to 
the high bending energy of stiff polymers. The increase in 
asphericity is continued in the semiflexible regime until it 
starts to decrease slowly for flexible polymers, clarifying 
ambiguous simulations [24, 25 . For phantom polymers 
the asphericity decreases down to (A) = 0.2625 in the 
Gaussian limit [26]. The three regimes of the phantom 
polymer resemble results for three-dimensional phantom 
rings, which, however, show a more pronounced decrease 
in asphericity in the flexible regime as the third spatial 
dimension enables more compact configurations. 

Focussing on the magnitude of the asphericity in the 
two dimensional models of phantom polymers and poly- 
mers with finite thickness, we find that the curves com- 
mence to deviate from each other in the semiflexible 
regime, showing fundamentally different values of as- 
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FIG. 4: The asphericities of the two models deviate in the 
semiflexible regime due to the influence of finite thickness. 
Through enhanced correlations between adjacent segments 
and prominent growth of the minor principal axis forced by 
pushing overlapping segments apart finite thickness makes the 
elliptical shape of ring polymers rounder. 



phericity in the flexible region. Hence, the ratio of the 
principal axes is strikingly different and the growth in size 
induced by the polymer's thickness is anisotropic along 
the principal axes yielding more spherical configurations 
compared to a phantom polymer. This is opposite to 
simulation results for open random walks in three di- 
mensions, where self- avoidance has been found to lead to 
slightly more aspherical configurations [27l [28] . In fact 
in three dimensions, a random walk is rarely intersect- 
ing its own trajectory. Then our present work suggests 
the remaining governing impact of self-avoidance in three 
dimensions to be an effective stiffening. Such a stiffen- 
ing yields more aspherical shapes for open semiflexible 
polymers in three dimensions. 

Confinement, like in the present case of two dimen- 
sional ring polymers, on the other hand, forces polymer 
segments to overlap much more frequently. Concerning 
two dimensional polymer rings, the notion of an aspheri- 
cal shape indicates that one principal axis is much longer 
than the other like in an ellipse. In the apices of the el- 
lipse, the segments are prone to overlap with neighboring 
segments on a local level, while segments in the convex 
part of the ellipse tend to overlap with segment being 
separated approximately half the total length along the 
contour. The finite thickness now effectively stiffens the 
polymer inducing less bending at the apices and it in- 
creases the minor principal axis by pushing segments in 
the convex region apart, see Fig. [I] This results in a 
more spherical configuration for polymer rings with fi- 
nite thickness, as observed in Fig. [4j As the asphericity 
distribution is highly skewed our statistical errors un- 
derestimate the actual error justifying slight deviations 
between our sets of data. Considering two dimensional 



polymer rings the good agreement between simulations 
and experimental data over the broad range of flexibili- 
ties manifests the role of finite thickness and its effects 
of effective stiffening and anisotropic change in shape. 

In conclusion, we have analyzed the impact of excluded 
volume caused by the finite thickness of polymer rings in 
two dimensions over the whole range of flexibility, both 
by observing DNA rings on mica surface and by com- 
puter simulations of phantom and finite thickness poly- 
mers. We find that the experimental data can only be 
explained by the latter, where each segment of the poly- 
mer is represented by an impenetrable tube. From the 
comparison of the different models, we determine two 
effects of finite thickness. Firstly, tangent-tangent cor- 
relation shows an enhanced correlation due to the steric 
constraints of the neighboring segments, leading to an 
effective stiffening observed in the semiflexible regime of 
the mean-square radius of gyration. Secondly, in the 
flexible regime Flory's swelling is recovered. However, 
the asphericity discloses an anisotropic change in shape. 
Manifesting these properties should enable a new under- 
standing of the conformation statistics of biopolymers 
such as DNA and F-actin. A basis on which biopolymer 
assemblies can be designed to develop new nanomateri- 
als. 

F.D., K.A. and E.F. acknowledge financial support 
of the German Excellence Initiative via the program 
"Nanosystems Initiative Munich (NIM)", of the DFG 
through SFB 486, and of the LMUinnovativ project 
"Functional Nanosystems (FuNS)". K.A. acknowledges 
funding by the Studienstiftung des deutschen Volkes. 
G.W. and G.D. acknowledge the support by the Swiss 
National Science Foundation through grants Nr. 200020- 
116515 and 200020-125159. 



[1] Seeman, N. C. Nature 2003, 421, 427-31. 
[2] Chen, J.; Seeman, N. Nature 1991, 350, 631-633. 
[3] Bloomfield, V. A. Biopolymers 1991, 31, 1471-81. 
[4] Hud, N. V.; Downing, K. H. P Natl Acad Sci USA 2001, 
98, 14925-30. 

[5] Lodish, H. F.; Berk, A.; Kaiser, C. A. Molecular cell bi- 
ology; Palgrave Macmillan: London, 2008; p 973. 

[6] Tang, J.; Kas, J.; Shah, J.; Janmey, P. Eur Biophys J 
2001, 30, 477-484. 

[7] Limozin, L.; Sackmann, E. Phys Rev Lett 2002, 89, 
168103. 

[8] Claessens, M. M. A. E.; Bathe, M.; Frey, E.; 

Bausch, A. R. Nat Mater 2006, 5, 748-753. 
[9] Cohen, A.; Mahadevan, L. P Natl Acad Sci USA 2003, 

100, 12141-12146. 
[10] Lau, A. W. C; Prasad, A.; Dogic, Z. Europhys Lett 2009, 

87, 48006. 

[11] Alim, K.; Frey, E. Phys. Rev. Lett. 2007, 99, 198102. 
[12] Flory, P. J. Principles of polymer chemistry, Cornell 

Univ. Press: Ithaca, N.Y., 1953; p 672. 
[13] Rubinstein, M.; Colby, R. H. Polymer Physics; Oxford 



5 



Univ. Press: Oxford, 2003; p 440. 
[14] Rivetti, C; Guthold, M.; Bustamante, C. J. Mol. Biol. 

1996, 264, 919-32. 
[15] Witz, G.; Rechendorff, K.; Adamcik, J.; Dietler, G. Phys. 

Rev. Lett. 2008, 101, 4. 
[16] Maier, B.; Radler, J. Phys. Rev. Lett. 1999, 82, 1911- 

1914. 

[17] Valle, F.; Favre, M.; De Los Rios, P.; Rosa, A.; Dietler, G. 

Phys. Rev. Lett. 2005, 95, 158105. 
[18] Ercolini, E.; Valle, F.; Adamcik, J.; Witz, G.; Metzler, R.; 

De Los Rios, P.; Roca, J.; Dietler, G. Phys. Rev. Lett. 

2007, 98, 058102. 
[19] DesCloizeaux, J.; Jannink, G. Polymers in solution', 

Clarendon Press: Oxford, 1990; p 896. 
[20] Aronovitz, J.; Nelson, D. J. Phys. (France) 1986, 47, 



1445-1456. 

[21] Stigter, D. Biopolymers 1977, 16, 1435-1448. 

[22] Rybenkov, V.; Cozzarelli, N.; Vologodskii, A. P. Natl. 

Acad. Sci. U.S.A. 1993, 90, 5307-5311. 
[23] Alim, K.; Frey, E. Eur. Phys. J. E 2007, 24, 185-191. 
[24] Camacho, C.; Fisher, M.; Singh, R. J. Chem. Phys. 1991, 

94, 5693-5700. 

[25] Norman, R.; Barker, G.; Sluckin, T. J. Phys. II 1992, 2, 
1363-1376. 

[26] Diehl, H.; Eisenriegler, E. J. Phys. A 1989, 22, L87-L91. 
[27] Cannon, J.; Aronovitz, J.; Goldbart, P. J. Phys. I 1991, 
1, 629-645. 

[28] Bishop, M.; Saltiel, C. J. Chem. Phys. 1988, 88, 3976. 



